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1. Introduction 

There are numerous distinctive features that characterize the behaviour of many-body 
systems with long-range interactions, features that are not present in systems with 
short-range interactions. These peculiarities concern both the equilibrium properties, 
such as the inequivalence of ensembles and negative specific heats in the microcanonical 
ensemble, and the out-of-equilibrium dynamical behaviour, such as the existence of long- 
lived quasi-stationary states and of out-of-equilibrium phase transitions. The study of 
these properties is interesting in its own, but it is also justified by the many different 
physical systems in which long-range interactions play the prominent role, e.g., self- 
gravitating systems [TJ |2j |3], unscreened Coulomb systems [I], some models in plasma 
physics [5] and in hydrodynamics [6] , and trapped charged particles [7] . Recent reviews 
give the state of the art of the subject [El [9]. 

In this paper, we treat the subject of the long-lived quasi-stationary states 
(QSS). They are out-of-equilibrium states in which the distribution functions are non- 
Boltzmannian, and their lifetime increases with the size of the system as given by the 
number N of degrees of freedom; this increase generally scales as a power law in N but 
it can also be exponential. It has to be emphasized that the QSS's are not related to 
the usual metastable states that are found also in short-range systems. The latter are 
realized by local extrema of thermodynamical potential (e.g., they are local maxima of 
the entropy or local minima of the free energy, if these quantities are computed as a 
function of an order parameter of the system), in which the system is trapped until it is 
driven away by some perturbation, and then heads towards the global extremum, i.e., the 
equilibrium state. Global and local extrema, i.e., equilibrium and metastable states, are 
obtained on the basis of the usual Boltzmann-Gibbs statistics. The evaluation of these 
states can be done following different routes; e.g., one can compute the partition function 
of the iV-body system or work at the level of the one-particle distribution function, but 
only with a direct relation to the Boltzmann-Gibbs statistics, that governs equilibrium. 

On the other hand, QSS's in principle have nothing to do with some sort of 
equilibrium state, global or local, of the system; nevertheless the system can be trapped 
for macroscopic times in these states. We want to underline that this fact does not 
imply at all any failure of the Boltzmann-Gibbs description of the equilibrium states. It 
simply means that the approach to equilibrium, that both in long-range and short-range 
systems should be described, at least approximately, by a kinetic equation, happens often 
in a manner, when long-range interactions are present, in which the system resides for 
macroscopic times in dynamical states that are very far from the equilibrium states. 

From this it should be clear that the ultimate reason for the existence of the QSS's 
should be looked into dynamics, i.e., into the representation of the dynamics via a 
kinetic equation. It turns out that the dynamics of many-body long-range systems can 
be described with a very good approximation, in a certain time range, by the Vlasov 
equation for the one-particle distribution function. This equation is also sometimes 
called the collisionless Boltzmann equation, since it represents the interactions between 
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the particles (whatever they might be, they can even be stars or galaxies in astrophysical 
problems) by a mean- field term, i.e., the two-body interaction potential averaged over 
the whole system. The theoretical justification for this fact can be looked upon at 
different levels of mathematical rigour. It is not our task here to give this justification, 
but we think it is useful to give a flavour of the reason from a physical point of view. 

If a system is initially prepared in a state away from equilibrium, it will evolve 
towards the equilibrium state because of the interactions between the particles. In 
short-range systems each particle interacts only with the nearby particles, and therefore 
the dynamical evolution is determined by the "collisions" of any given particle with the 
few others surrounding it. Since the position of nearby particles are strongly correlated, 
it is completely useless to approximate the field acting on a particle by an averaged field, 
but it is necessary to find a way to describe, in a kinetic equation, the collisions between 
close particles. In long-range systems, even in the cases where the field at close distances 
is strong, the field acting on a particular particle is determined by all the others, and it 
looks quite plausible that an averaged field can be a good approximation. Obviously, at 
the end, the "collisional" regime will take its toll and the Boltzmann-Gibbs equilibrium 
will be realized, but before that point, the dynamics, to a high degree of approximation, 
will follow the evolution determined by the Vlasov equation. 

It is then not surprising that a considerable amount of work has been dedicated 
to the properties of the Vlasov equation in relation with important systems like self- 
gravitating and plasma systems. It is worth noting a difference between these two 
classes of systems. In fact, globally neutral plasmas have a spatially homogeneous 
equilibium state, and it is in reference to this case that the theory of the Vlasov 
linear stability of small perturbations of the homogeneous equilibrium state has been 
developed [TO] [TT] [12] . On the contrary, the necessarily inhomogeneous states of large 
but finite self-gravitating systems has motivated the research on the stability properties 
of inhomogeneous stationary states of the Vlasov equation [U HH [151 EE ITT] . 

The Hamiltonian Mean Field (HMF) model [T8l [19] . a simple ID toy model 
of systems with long-range interactions, has been very useful to study the various 
statistical and dynamical properties of long-range systems. Also, a good amount 
of work has been dedicated recently to the stability of Vlasov stationary states. 
The comments of the previous paragraph suggest that one should be interested in 
the stability of both homogeneous and inhomogeneous states. However, mainly 
homogeneous states have been considered (both theoretically and numerically, see, e.g., 
Refs. [191 EHl HU [221 [231 I2H I25J [26]), and there exists only few theoretical results 
for inhomogeneous states (see Sees. 3.3 and 4.4 of [24], Appendix F of [26], and 
[27] for polytropes). In this paper, we present a criterion for the Vlasov stability of 
inhomogeneous stationary states of the HMF model that generalizes the known criterion 
valid for homogeneous states. The results presented are based on the stability conditions 
for Vlasov stationary states that have been derived in the (essentially astrophysical) 
literature. We rederive here the results, and extend them, leaving the interaction 
potential unspecified. We will introduce the HMF interaction potential when the 
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mentioned results will be used to obtain explicit necessary and sufficient conditions 
on the stationary distribution function. 

We emphasize that the stability conditions already in the literature are in the form 
of relations that have to be satisfied by the perturbations to the stationary states. To 
our knowledge, such conditions have not been yet transformed in explicit conditions 
that have to be obeyed by the stationary distribution itself. Therefore, such explicit 
conditions are the core of this paper. 

It turns out that the necessary and sufficient conditions of linear stability can be 
simplified at the price of obtaining conditions that are only sufficient. This can be done 
in the framework of the formal stability of the Vlasov stationary states. Therefore, 
for more completeness and for a useful comparison, we find it useful to treat also the 
problem of formal stability. 

Section [2] introduces the stationary states in the framework of the Vlasov equation. 
Sections from [3] to [6] derive results later employed for the study of linear stability and 
formal stability. This is done, for the HMF model, in Sections [7] and El Section [9] 
contains the conclusions. 



2. The Vlasov equation and the quasi-stationary states 

This paper will use the HMF model as a benchmark for our analysis, and therefore it 
is convenient to introduce the Vlasov equation from the beginning in this framework. 
Consider N particles of unit mass moving on a circle, with the Hamiltonian of the system 
given by: 

a) 

i=i i^j 

where Q\ G [0, 2tt] is the angle giving the position of a particle on the circle and 
— oo < Pi < oo is its linear momentum (equal to the velocity since the mass is 
unitary). The normalization of the interaction potential V(9i — 9j) is the usual 
one introduced in order to have an extensive energy; it is equivalent to a system-size- 
dependent rescaling of time, and it does not affect the study of the properties of the 
system. The Vlasov equation associated to this system, governing the evolution of the 
one-particle distribution function f(9,p,t) is: 

df(e, P ,t) df(o,p,t) d$(o,t;f)df(e, P ,t) _ n 

dt +P 39 89 dp - ' [) 

where $(<9,i; /) is the mean field potential: 

/oo p2tt 
dp' / d9'V(9-9')f(9',p',t). (3) 
-oo JO 

The last equation shows that the Vlasov equation (j2j) is a nonlinear integrodifferential 
equation. It is immediate to see that it conserves the normalization of f(9,p,t): 

dp d9f(9,p,t) = l, (4) 



Stability of inhomogeneous quasi- stationary states 



5 



and the total energy, given by: 

e= r dp C de + ^ {d,t] f) ) f{e,p,t) ■ (s) 

The HMF model is obtained when V(9 — 9') is a cosine potential. In the following, we 
will present the criteria that determine the dynamical stability of stationary states of 
the Vlasov equation (J2J) for a general potential V(9 — 9'). These criteria can be trivially 
generalized to arbitrary space dimension. Afterwards, we will specialize to the HMF 
model when the relations will be transformed in explicit conditions for the stationary 
states. 

As we have described in the Introduction, there is a time regime in which the 
dynamics of the iV-body system is, to a high degree of accuracy, represented by the time 
evolution of the one-particle distribution function as governed by the Vlasov equation. 
In this framework, it is natural to expect that the QSS's will be associated to stationary 
states (i.e., time independent states) of this equation. Of course, these stationary states 
should also be stable, i.e., a small perturbation should not drive the system away from 
the stationary state. Then, it is natural to be interested in the dynamical stability of 
the stationary states of the Vlasov equation. It is easy to see that any function of the 

2 

form f(9,p) = /(y + $(6*; /)) is stationary; therefore, in principle one is interested in 
determining the stability of any such function. 

2 

The function /(y + <I>(0; /)) is by definition linearly stable if it is possible to choose 
the norm of the perturbation 5f(9, p, t) at time t — such that this norm remains smaller 
than any (small) positive number, provided that the dynamics of the perturbation is 
governed by the linearized Vlasov equation, with the linearization made around the 
stationary state. Introducing the individual energy: 

e(0,p) = £ + mf)> (6) 
the linearized Vlasov equation for 5f(9,p,t) is easily obtained as: 

d5f(9,p,t) dSf(9,p,t) d$(fl;/) dSf(9,p,t) d$(9,t;6f) 

dt + p —9 d9 dp d9— pf = ° ' (?) 

where the potential $(#; /) is constant in time (and thus the partial derivative with 
respect to 9 has become a total derivative) and where in the last term the functional 
dependence of / on (9,p) only through e(9,p) has been exploited. The problem of 
the linear stability associated to this equation has been treated long ago by Antonov 
[T3] in astrophysics. In order to have a self-contained presentation, in section H] we 
will reproduce, although in a somewhat different formulation and scope, the Antonov 
results, that later will be used to derive the stability conditions on the stationary state 
/(y + 3>(0; /))• These will be necessary and sufficient conditions of linear stability. 

As we have previously underlined, we are also interested in less refined stability 
criteria, that provide only sufficient, but simpler, conditions of linear stability. To that 
purpose, we can use the notion of formal stability of a stationary point of a general 
dynamical system. The stationary point is said to be formally stable [28] if a conserved 



Stability of inhomogeneous quasi- stationary states 



6 



quantity has an extremum at the stationary point and if the second variation about 
the stationary point is either positive definite or negative definite. It can be shown [28] 
that formal stability implies linear stability; therefore, proving that a stationary point 
is formally stable gives a sufficient condition for linear stability. Formal stability is also 
a pre-requisit for nonlinear stability, although formal stability does not imply nonlinear 
stability for infinite dimensional systems. With a slight extension of definition, we will 
consider the formal stability of a stationary point also for the cases in which both the 
extremization of the conserved quantity and the sign of its second variation are studied 
under some constraints. The problem will be related, as it will be clear, by the fact 
that if one finds that the second variation has a definite sign for the unconstrained case, 
this is sufficient to have the same definite sign also for the constrained case. However, 
the converse is wrong and this is similar to the notion of ensembles inequivalence in 
statistical physics [29J. 

To study the formal stability problem, we use a general result obtained in Ref . [30] . 
Again for a self-contained presentation, we find it useful to briefly reproduce it here. 
This is done in the following section, in slighlty more general terms than those necessary 
for the linear Vlasov equation. 



3. The maximization of a class of functionals 

Let £ be the functional space of real differentiable functions f(6,p) defined for < 9 < 
2it and p £ 1Z. We also assume that the functions decay sufficiently fast for p — » oo so 
that the integral of p 2 f is finite. The scalar product in this space is naturally defined 
by: 

/OO /*27T 
dp / de g (e,p)f(e, P ). (8) 
-oo JO 



We consider here the problem [30] of finding the constrained maximum of the functional: 

/OO /*2"7r 
dp / d0C(/(M) , (9) 
-oo Jo 

with the function C(x) at least twice differentiable and strictly convex, i.e., with the 
second derivative strictly positive. The constraints are given by two functionals. The 
first is a linear-quadratic expression: 

/*oo r2"K 

E= dp dea(6,p)f(9,p) 



oo JO 

oo r-2-K /-oo r2n 



-r /*oo pzn poo pzir 

+ - dp d9 dp' d9'f(9,p)b(9,p,9',p')f(9',p'), (10) 

^ J-oo Jo J-oo Jo 

that we can call the "total energy" (in analogy with Eq. fl5])), and that is constrained 
to have a given value Eq. In this expression, a and b are two given functions, with b 
possessing the symmetry property b(9',p',9,p) = b(9,p,9' ,p'). The second constraint is 
the normalization: 

poo r2ir 

1= dp d9f(9,p) = l. (11) 
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By using Lagrange multipliers (3 and fi the extremum is given by equating to zero the 
first order variation: 



5S - /35E - fi5I = 

We thus have: 

-C(f(e, P ))-p 



2tt 



a(6, P )+ / dp' / de'b(e,p,e',p')f(e', P ') 



- fl = 



If we denote the "individual energy" (in analogy with Eq. (Q) by: 



2- 



(12) 
(13) 

(14) 



(15) 



e(9,p) = a(9,p) + / dp' / d9' b(9,p,9',p')f(9',p') , 

J ~oo JO 

then the extremum relation can be written as: 

-C'(f(e,p)) = 0e(9,p)+ti. 
From the convexity property of C(x) it follows that this relation can be inverted to give: 

f(9,p) = F(J3e + »)=f{e), (16) 

where F is the inverse function of —C. Inserting this function in Eqs. f TTOT) and (11 1H 
we obtain the values of the Lagrange multipliers. It is clear that Eq. (1141) is also 
a consistency equation. We note that in order to interpret f(0,p) as a distribution 
function, acceptable functions C(f) in Eq. ([9]) are only those that, through Eq. (fl6|) . 
provide a positive definite function. We also note the identity: 



f(e(9,p)) 







(17) 



C"(f(6,p)Y 

that is obtained by differentiating Eq. (TIB]) . From the convexity property of C(x) 
it follows that ^f'(e(9,p)) is negative definite. On the other hand, since we have 
-jj- = pf'(e(9,p)), the integrability in p requires that, if f'(e(9,p)) has a definite sign, 
this sign must be negative. Therefore, (3 is restricted to positive values. We therefore 
conclude that the extremization of S at fixed E and / determines distribution functions 
of the form / = f(e) with f'(e) < 0. 

The extremum so obtained will be a maximum only if the second order variation 
of S[f], for all the allowed displacements 8f(9,p), is negative definite. The variation of 
the functional S[f] is given, up to second order, by: 

f27T 



6S 



dp 



d9 



o 



C'(f(9,p))8f(9,p) + -C"(f(9,p))(8f(9,p)) 2 



(18) 



with the derivatives computed at the extremal point. Then, from Eqs. ( I15p . (I16p and 
f fT7|) we obtain: 

-2tt 



5S 



dp 



d0 



2f'(s(9,p)) 



(19) 
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We can transform this expression by using the fact that the variations of E and I must 
identically vanish for the allowed displacements 5f(8,p). These variations are given by: 

POO p2lT 

5E= I dp I d9e(9,p)5f(9,p) 

+ 
and 



oo 

-1 POO p2n POO p2n 

-/ dp de dp' de , 5f(e, P )b(e,p,9',p')5f(e , ,p , ) = o (20) 

^ J-00 JO J-00 Jo 

/oo p2ir 
dp / d65f(6,p) = (21) 
■00 Jo 

respectively. Adding to Eq. ( fl9l) the zero valued expression —j35E — fi5I we have: 

poo p2tt 



-i roc r27T roo r2n 

--pi dp d9 dp' d6'5f(6,p)b(8,p,6',p')5f(6',p'). (22) 

^ J-oo JO J -oo JO 

The right-hand side of this expression must be negative definite for all allowed 
displacements 5f(9,p), i.e., for all those that at first order do not change E and /. 

The problem of maximizing S[f] at constant E and I can be shown to be equivalent 
to that of minimizing the energy E at constant S and I. In fact, using Lagrange 
multipliers 1/(5 and —fi/(3, the equation of the first order variation is now: 

5E - + ^51 = 0, (23) 

which is the same as Eq. (fT2~|) ; then the solution is again given by Eq. (IT31 . Now, we 
have to study the variation of E, that up to second order is given by the left-hand side 
of Eq. (JUD, i.e.: 

dp / d6e(6,p)5f(6,p) 

-oo JO 



-oo 

-i roo /*27r ^oo /*2tt 

+ -/ dp d6 dp' d6'5f(6,p)b(6,p,6',p')5f(6',p'). (24) 

^ J-oo Jo J-oo Jo 

The variations of S and / must identically vanish; the latter is expressed by Eq. (121 p . 
while the former is given by: 



5S = / dp de 



2- 



(Pe(e,p) + fi)5f(9,p) + l l* (8f{9,p)? 



ee . (25) 



Adding to Eq. ( 123]) the zero valued expression —\o~S + %dl we arrive at: 

-i poo p2n -i 

5E=-- dp d6———(5f(8,p)f 
2 J-oo Jo / 

-| poo r2n poo p2n 

+ - dp d6 dp' d9'5f(9,p)b(9,p,6',p')5f(9',p'). (26) 

^ J-oo Jo J-oo Jo 

This is the same as Eq. (1221) divided by —ft < 0. Therefore, Eq. (|26j) is positive definite, 
i.e., £7 is minimum at the stationary state, if Eq. (1221) is negative definite. To complete 
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the proof of the equivalence between the maximization of S at constant E and I and 
the minimization of E at constant S and / we have to see that in both cases the allowed 
displacements 5f(8,p) are the same. This can be deduced in the following way. For 
Eq. (1221 the allowed displacements are all those that at first order give SE = 51 = 0. 
By Eq. (fl2l) . they also at first order give 5S = 0; then, they are also allowed for Eq. 
f )26|) . In turn, the allowed displacements for Eq. (126]) are all those that at first order 
give 5S = 51 = 0. By Eq. ([23]) . they also at first order give 5E = 0; then, they are also 
allowed for Eq. ( )22|) . This concludes the proof. 

After treating the problem of the linear stability of the stationary states of the 
Vlasov equation, we will use the results of this section to study their formal stability. 



4. The linear stability of Vlasov stationary states 



The energy functional in the previous section, Eq. ( flO|) . was more general than the one 
associated to the Vlasov equation, Eq. ([5]). The former reduces to the latter when the 
functions a(9,p) and b(9,p,9',p') are related to the kinetic energy and to the potential 

2 

energy of the system, respectively; namely, when a = y and b = V{6 — 9'). For 
convenience, we rewrite here the relevant expressions. We have the individual energy 

e(M = y + *(*;/) , (27) 
with the mean field potential 

/OO /*27T 
dp' / d9'v(9-9 , )f(9 , ,p'), (28) 
-oo JO 

and the total energy 

E[f] = j°° dp j\e 0j + h(6,t; f)j f($,p,t) . (29) 

We will consider stationary states that are associated to the extremization of functionals 
of the form ([9j). We have seen that in this case f(0,p) = f(e(9,p)) with f'(e(9,p)) < 0. 
For more compactness, we will use the notation / ~f(9,p) = f'(e(9,p)). As remarked 
above, we follow and complete the treatment of Antonov [13] . 

In the following we will need the usual extension of the scalar product defined in 
Eq. OH]) to complex valued functions: 

/oo r2TT 
dp / &0g*{0,p)f{0,p), (30) 
-oo ^0 

where the asterisk denotes complex conjugation. 

It is not difficult to see that the linearized Vlasov equation (J7J), determining the 
dynamics of 5f(9,p), can be cast in the form: 

^■r(0,P,t) = -i{9,p){DK5j){9,p,t) , (31) 
where D is the antisymmetric linear differential operator (advective operator): 

(Dg){9,p) = P^9(e,p) - ^^9(0, p) , (32) 
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while K is the linear integral operator: 

(Kg)(9,p) = -L-g(e,p) - m 9) ■ (33) 
7(0, p) 

This can be obtained by seeing that De = 0; then, the action of D on any function of 
e gives zero; in particular D'y = 0. Finally, it can be easily checked that the operator 
B = DKD, needed shortly, is hermitian. 

The stationary point of the dynamics Sf(9,p) = will be linearly stable iff: (i) all 
nonsecular solutions of the type 5f(9,p, t) = 5f(9,p, 0) exp(At) have eigenvalues A with 
non positive real part; (ii) in the presence of secular terms (i.e., if there are eigenvalues 
with an algebraic multiplicity larger than the geometric multiplicity), the associated 
eigenvalue must have a negative real part. 

Before proceeding further, we need to put in evidence the properties of the operators 
D and K when acting on functions that are either symmetric or antisymmetric in p. 
From the definition of D in Eq. ( 1321) it is clear that D transforms symmetric functions in 
antisymmetric functions, and viceversa. Concerning K, we see from its definition in Eq. 
( 1331) that it maintains the symmetry of the functions; however, since for antisymmetric 
functions g a (9,p) we have $(0; g a ) = 0, the action of K in this case simplifies in 

(Kg a )(9,p) = —}—g a (9,p). (34) 
7(0, p) 

We now suppose that 5f(9,p; A) is the eigenfunction associated to the eigenvalue 
A. We then have, from Eq. fl3"TT) : 

X5f(9,p;X) = - 1 (DK5f)(9,p;X), (35) 

where for simplicity we have dropped the dependence of 7 on the coordinates. We now 
separate 5f in the symmetric and antisymmetric parts: 5f = 5f s + 5f a . Taking into 
account the mentioned properties of the operators D and K we obtain 

X5f s (9,p;X) = -(D5f a )(9,p;X) (36) 

and 

X6f a (9,p; A) = - 7 (DK5f s )(9,p; A) . (37) 

If we multiply the second of these equations by A and substitute X5f s from the first we 
have: 

X 2 6f a (9, P] A) = ^(DKDSf a )(9,p; A) = -y(B6f a )(e,n A) . (38) 

If A ^ 0, Eqs. ( |36l) and ( 1371) prove two properties. The first is that an 
eigenfunction cannot be either symmetric or antisymmetric, but both components must 
be nonvanishing. The second is that, if 5f s + Sf a is associated to the eigenvalue A, then 
£/s — fifa is associated to the eigenvalue —A. 
From Eq. ( l38|) we have: 



i2 6fa(0,p;\) 



A 2 ' =(B5f a )(9,p;X). (39) 

7(0, P) 
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The scalar product of both sides of this expression with Sf a gives: 

A 2 f°dp [ 2 \e lSfa{9 ( f>^ 2 = (5f a ,B5f a ), (40) 

while the scalar product of its complex conjugate with 8f* gives, exploiting the 
hermiticity of B: 

(A*) 2 dp r d e l<5/a(g /\ A)|2 = itf/ fl > • (41) 



7 (0,p) 

Since 7 is negative definite, Eqs. (14"0~]) and (j4Tj) imply that, if A 7^ 0, A 2 is necessarily real, 
and therefore A is either real or pure imaginary. In much the same way, it can be shown 
that, if 8f(9,p;Xi) and 5f(9, p; A 2 ) correspond to two different non zero eigenvalues, 
then: 

(5f a (X 1 ),B5f a (X 2 )) = Q. (42) 

The case A = will be considered later. 

Let us now assume that an eigenvalue A has an algebraic multiplicity larger than its 
geometric multiplicity. Then, if 5f(9, p; A) is an eigenvector associated to this eigenvalue, 
there will also exist a solution of Eq. fl3T|) given by [tSf(6,p; A) + 5f^(9,p; A)] exp(Ai). 
Substituting in Eq. ( |3T|) and using Eq. ( 1351) we obtain: 

5f(6,p; A) + \5fU(9,p; A) = - 7 (DK5f^)(9,p; A) . (43) 

Separating both Sf and 5f^' in the symmetric and antisymmetric parts we obtain: 

8f s (9,p; A) + X5f^(9,p; A) = -(D5f^)(9,p; A) (44) 

and 

5f a (9,p; A) + X5fU(9,p; A) = ->y(DK6fP)(0,n A) . (45) 

If we multiply the second of these equations by A and we substitute Xd~fs^ from the first 
we have: 

X5f a (9,p; A) + X 2 5fi 1 \9,p; A) = j(B5f^)(9,p; A) + j(DKSf s )(9, p; A) . (46) 

Using Eq. ()3Tj) we arrive at: 

2Ai/.[», B »Utf(»,,A| <„ 
7(0, p) 

The scalar product of this expression with 8f a gives: 

- 2A|5/ a (0,p; A)| 2 + A 2 5/ a *(0,p; A)5/i 1} (0,p; A) 



(4* 



i-00 Jo 7(0, P) 

Substracting from this equation the one that is obtained by forming the scalar product 
of the complex conjugate of Eq. ( 1391 with 5 fa we have: 

A fdp r de \lM^l = . (49) 

J -co Jo 1\?,P) 
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If A 7^ we deduce that 8f a (9,p;X) = 0, since 7 is negative definite; then also 
Sf s (9,p;X) = and thus 5f(9,p;\) = 0. This shows that no eigenvalue A different 
from has an algebraic multiplicity larger than its geometric multiplicity. 

We now consider the case A = 0. We note that the presence of a zero eigenvalue 
with an algebraic multiplicity larger than the geometrical multiplicity would imply 
the presence of a solution of the equation of motion (|3T|) of the form [td~f(9,p;0) + 
Sf^(9,p; 0)], and therefore the stationary point 5f(9,p) = would be linearly unstable. 
In the following we need to consider also the possibility that the difference between the 
algebraic and the geometric multiplicity of the zero eigenvalue is such that solutions 
with higher powers of the time t exist. 

For A = 0, Eqs. fl36|) and fl37j) become: 

(D5f a )(9,p;0) = (50) 

and 

(DK6f s )(9,p;0) = 0. (51) 

We note that Eq. ([50]) implies {B5f a ){9,p; 0) = 0. 

In the case the eigenvalue A = has an algebraic multiplicity larger than the 
geometric multiplicity, we obtain, from Eqs. (1441) and (T4"5]) : 

Sf s (9,p;0) = -(D5f^)(9,p;0) (52) 

and 

5f a (6,p;0) = - 7 (DK6fW)(9,p;0). (53) 

Substitution of Eq. ([52]) in Eq. (J5T) gives (B5f^ ] )(9,p;0) = 0. Dividing both sides of 
Eq. (]53l by 7 and forming the scalar product with 5f a we get: 

r dp f\ 6 ^ U0 ( f^ )? = -(6f a ,DK6fV) = (D6f a ,K6f?>) = 0, (54) 
J -00 Jo 7(0, P) 

where in the last step we have used Eq. ( 1501) . It follows that Sf a (9,p; 0) = 0, and then 
5f s (6,p; 0) 7^ in order to have a non trivial solution. Then, from Eq. f[52l) we obtain 
that 5jj 1) {9,p; 0) + and {D5fj l 1) ){9,p; 0) = -5f s {9,p; 0) + 0. 

If there exists a solution of the equation of motion ( 13T1) of the form [t 2 Sf(9,p; 0) + 
tSf ( - 1 \9,p; 0) + 6fW{6,p; 0)], then the evaluation of (I3T1) at t = gives: 

5f^(9,p;0) = -y(DK5f®)(e,p;0). (55) 

The usual separation in the symmetric and antisymmetric parts gives: 

6f!; 1 \9 1 p;0) = -(D5fW)(9 1 p;0) (56) 

and 

SfPfrPiO) = -l(DK6f^)(9,p;0). (57) 

The substitution of Eq. ( [56]) in Eq. ( [52]) . taking into account that Sf a (9,p; 0) = 0, 
gives (B6f! 1 2) )(9,p;0) = 0. Finally, it can be shown that a solution of the form 
[t 3 5f(9,p; 0) + t 2 5f^(9,p; 0) + t6fW(6,p; 0) + 6fW(6,p; 0)] cannot exist. In fact, the 
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evaluation of ( 13 ip at t = gives, after separation in the symmetric and antisymmetric 
parts: 

6fP(6,p ] 0) = -(D5f^)(6,p;0). (58) 
Substitution in Eq. (jSTj) gives: 

5/«(^p;0) = 7 (^/i 3) )(^p;0). (59) 

Dividing both sides by 7 and forming the scalar product with 5 fa we obtain: 

r dp f 2 \ \6fP{0,pM> = {Sf (i) )BSf p } = { B5fU,8fW) = 0, (60) 
J-00 Jo 1\?,P) 

since (B8fa)(9,p; 0) = 0. This implies that 5fa (9,p; 0) = 0, that in turns gives, from 
Eq. fl52|) . 5f s (9,p; 0) = 0. This is not acceptable, since we already have 5f a (9,p; 0) = 0. 

Summarizing all the results, we have that the initial value of 5f(8,p) can be 
decomposed in general as: 

8f(e,p)= ^c(\)5f(e,p-x) 



+ 



J2[cjSf J (e i p;0) + c i j L) 6f^\e,p;0) + cf5fP(e,p;0)] , (61) 
3=1 

where the sum over the nonzero eigenvalues stands also for an integral if the eigenvalues 
are continuously distributed. The sum over r different contributions coming from the 
zero eigenvalue takes into account possible separation into disjoint eigenspaces; some of 
the corresponding functions might have only 5fj(9,p;0) different from zero, and some 
only df^O) and 5 f^ (9 , p; 0) . The fact that (B6f j>a )(9,p;0) = (B8f^)(9,p;0) = 
(B5fj 2 J)(9,p;0) = implies that if we take the antisymmetric part of d~f(9,p) in Eq. 
fl6T|) and form the scalar product (6f a ,B5f a ), only the eigenfunctions corresponding 
to the eigenvalues A different from zero contribute. Precisely, using the orthogonality 
property (pf2l) we have: 

(8 fa, BSf a ) = \c(X)\ 2 (5f a (X), B5f a (X)} . (62) 
With the use of Eq. (HI]) we get: 



/OO PZ7T 
dp / d# 
-00 JO 



^ \5f a (9,p;\)\* 



1(9, p) 



(63) 



A^O 

This is the first important expression of this section. From it we deduce that, if the 
stationary point 8f(9,p) = is linearly stable, then necessarily the last expression is 
positive for all cases in which not all the coefficients c(A) are zero. In fact, the linear 
stability requires that all eigenvalues different from zero are pure imaginary (we recall 
that if a real negative eigenvalue exists, also its opposite exists and leads to instability; 
we also recall that the antisymmetric part of any eigenfunction corresponding to an 
eigenvalue different from zero is nonvanishing) . From this and from the fact that 7 is 
negative definite, our statement follows. 
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The scalar product (163 p can be zero for a nonvanishing Sf if one or some of the 
coefficients Cj, c*p are non zero. We know that, if one of the eigenspaces, e.g. the 
one corresponding to j = jo in Eq. ( ISTI) . has an algebraic multiplicity larger than its 
geometric multiplicity, the stationary point 5f(8,p) = is linearly unstable. In this 
case, taking 5f(6,p) = 5f^\9,p) we have (5f a ,B5f a ) = with, according to what was 
proved just after Eq. fl54|) . (DSf a )(9,p) ^ 0. On the other hand, if all the eigenspaces 
corresponding to A = have equal algebraic and geometric multiplicities, then only the 
terms with Cj appear in Eq. ( IBTj) as the contribution from the zero eigenvalue. But in 
this case, it follows from Eq. ( |50j) that the scalar product (163]) can be zero only for a 
function Sf(9,p) such that (D5f a )(9,p) = 0. 

In conclusion, the necessary and sufficient condition for the linear stability is that 
the scalar product ( |63l) is nonnegative, and it is zero only for functions Sf(9,p) such 
that (D5f a )(9,p) = 0. 

We now rewrite the scalar product (Sf a , BSf a ) in another way. Precisely: 

(5f a ,B6f a ) = (6f a ,DKD6f a ) = -(DSf a ,KDSf a ) 

= - (D5f a , -DSf a ) + (D5f a , $(D6f a )) , (64) 

7 

where use has been made of the definition of the operator K in Eq. f[3"3l . Using the 
definition of $ in Eq. ([3]) we finally obtain: 

/oo /*2"7r -i 

d P / de——\(D8f a )(6,p)\ 2 (65) 
Jo lK e iP) 

/oo p2tv poo flit 

dp d9 dp' / d9' (D5f a )*(9,p)V(9 - 9')(D5f a )(9',p') . 
-oo JO J -oo JO 

We should note two things. Firstly, since our equation of motion is real, it is always 
possible to choose Sf(9,p) real. Secondly, it is not difficult to check, form the last 
expression, that if we consider the whole function d~f(9,p) instead of its antisymmetric 
part 5f a (9,p), we will always have (5f,B5f) > (5f a ,B5f a ). We then arrive at the 
following necessary and sufficient condition of linear stability: 



dp \ d9——((D5f)(9,p)Y (66) 
7(0, p) 



/OO pZTT POO pZTT 

dp d9 dp' / d9'{D5f){9,p)V{9-9'){D8f){9',p')>0 
-oo ^0 J -oo Jo 

for any 5f(9,p), with the equality holding only when (DSf)(9,p) = 0. This is the main 
expression of this section. In treating the particular case of the HMF model, it will be 
the basis to obtain a condition on the stationary state f(e(9,p)). 



5. The energy principle and the most refined formal stability criterion 

The stability of stationary states of the Vlasov equation has been studied by Kandrup 
[T6| [T7] , introducing a Hamiltonian formulation for this equation, obtaining a sufficient 
condition of linear stability. We give here few details, mainly to show that, for stationary 
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states of the form f(s(9,p)), with f'(e(9,p)) < 0, this condition becomes identical to 
the one given in Eq. ( 166]) . thus becoming also necessary. 

This approach is based on a Hamiltonian formulation of the Vlasov equation, and 
on the observation that the Vlasov dynamics admits an infinite number of conserved 
quantities, called Casimir invariants, given by: 

/*OG /*27T 

C A [f]= / dp / d9A(f(9,p)), (67) 



J-oo JO 

for any function A(x). Note that functionals of the form ([9]) are particular Casimirs. 
The Hamiltonian formulation of the Vlasov equation ([2]) for f(9,p) is realized by casting 
it in the form 

% + iM = Q > (68) 



where e is the individual energy fl27|) and the curly brackets denote the Poisson bracket: 

da db da db 

M = wWp-Wpd§- (69) 

It can be shown that, if f{6,p) is a stationary state of the Vlasov equation of a general 
form, then a sufficient condition for its linear stability is the following [17] : the difference 
between the total energy (129]) computed at the perturbed state f(9,p) + Sf(9,p) and 
the total energy computed at the stationary state f(9,p) is positive for all Sf(9,p) 
that conserve all the Casimirs. In other words, f{9,p) is linearly stable if it is a local 
minimum of energy with respect to perturbations that conserve all the Casimirs. This 
forms the most refined formal stability criterion. These so-called "phase-preserving" or 
symplectic perturbations can be expressed in the form: 

f(9,p) + 8f(9,p) = eMf(9,p), (70) 

for some "small" generating function a(9,p). They amount to a re-arrangement of phase 
levels by a mere advection in phase space. Expanding to second order in a we have 

f(9,p) + 5f(9,p) = f(9,p) + {a, /} + l - {a, {a, /}} . (71) 

The corresponding expansion of the total energy f|29|) is easily obtained as: 

dp / d9e(9,p){aj}(9,p) 

-oo JO 



-i /*GO pin 

+ - dp d9e(9,p){a,{a,f}}(9,p) 

'-oo JO 
foo /•2tt poo /•2ir 



2 



-I pOO P Z7T pOO pZTT 

+ -/ dp d9 dp' d9' {a,f}(9,p)V(9-9'){a,f}(9',p'), (72) 

^ J-oo JO J-oo JO 

where the individual energy ( 127]) at the stationary distribution has been used (and for 
clarity the explicit dependence of the Poisson brackets has been written). Now, it is 
possible to exploit the identity: 

/oo p2iT poo r2ir 

dp / d^ Cl {c 2 , C3 }= / dp / d9c 2 {c 3 , Cl } (73) 
-oo JO J-oo JO 
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for any three functions Ci, 02,03. Using in addition that for a stationary state 
{/, e} = 0, we have that the first line of the right-hand side of Eq. f J72|) . i.e., the 
first order variation of the total energy, vanishes. This shows that the total energy 
at a stationary distribution is an extremum (5E = 0) with respect to symplectic 
perturbations. However, this does not guarantee that it is an extremum with respect to 
all perturbations. Using again the identity f!73|) . the second order variations of energy 
deduced from Eq. f!72|) are: 

-r poo p2tt 

5 2 E=- dp d9 {e,a}(6,p){a,f}(6,p) 



2 j-00 jo 

-1 f'OG p27V pOO f'l'K 

+ -/ dp de dp' de' {aj}{e,p)v{e-e'){a,f}{e',p'). (74) 

^ J-oo JO J-oo JO 

The positive definiteness of this expression is a sufficient condition of linear stability. 

We now suppose that the stationary distribution function is a function of e(9,p). 
In this case {a, /} = f'(e) {a, e}. Furthermore, we also have {a,e} = Da, where D is 
the linear differential operator defined in Eq. ( 13"2"j) . Using that De = 0, and therefore 
Df'(e) = 0, Eq. (!74|) becomes in this case: 



-| poo p2n poo p2-k 

- dp d9 dp d6' (Da){6,p)V{6-e'){Da){6',p') 

2 J-oo JO J-oo JO 



s2E -2L dr 'l d0 7WW) {{mmp>r (75) 

POO f2lT poo r2ir 

2 

where a = f'(e)a. The positive definiteness of this expression is exactly the necessary 
and sufficient condition of linear stability (1661) . recalling the definition of / ~f(6,p). 

It is convenient to introduce the notation Sf(6,p) = Da(6,p). Then, the necessary 
and sufficient condition of linear stability can be written 

''OO p2n -i 

dp / d9——(6f(6,p)) 2 (76) 

Jo n/{o,p) 

/oo / > 2tt poo p2ty 

dp de dp' de' 8f{e,p)v{e-e')8f{e',p') > o, 
oo Jo J~oo Jo 

for any perturbation of the form Sf(e,p) = Da(e,p) where a(e,p) is any function. These 
perturbations correspond to a mere displacement (by the advective operator D) of the 
phase levels, i.e. to dynamically accessible perturbations. It is straightforward to check 
by a direct calculation that these perturbations conserve energy and all the Casimirs 
at first order (this is of course obvious for symplectic perturbations). Indeed, using 
Sf(e,p) = {a, /} and identity (1731) . we get 

/oo p2n poo p2tt 

dp / d65f(e,p)e= dp de{a,f}e 
-oo Jo J-oo Jo 

/oo p2tt 
dp / d#{/,e}a = 0, (77) 
oo Jo 



and 



/oo p2it poo p2it 

dp d6 5f(9,p)A'(f)= dp de{a,f}A'{f) 
-oo Jo J-oo Jo 
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— oo 



■oo 



dp 




d9{f,A'(f)}a = 0. 



(78) 



6. Less refined formal stability criteria: sufficient conditions of stability 

As we have seen above, the minimization of E with respect to symplectic perturbations, 
i.e. dynamically accessible perturbations that conserve all the Casimirs, is a necessary 
and sufficient condition of linear stability. It is also the most refined criterion of 
formal stability since all the constraints of the Vlasov equation are taken into account 
individually. Less refined formal stability criteria, that provide only sufficient (albeit 
simpler) conditions of linear stability can be obtained by relaxing some constraints. 

6.1. The "microcanonical" formal stability 

As we have seen in section El a stationary state of the form f(0,p) = f(s(9,p)) with 
f'(e) < is obtained by extremizing a functional S of the form (Q, with the function 
C(x) given by Eq. (fT6l) . under the constraints given by the normalization /, Eq. ( fill , 
and the total energy E, Eq. ( 1291) : but also by extremizing the total energy E at constant 
S and /. Furthermore, we have proven that a maximum of S at fixed E and / is a 
minimum of E at fixed S and / and viceversa. Since we have seen that for distributions 
of the form f(e(6,p)) the minimization of E with respect to symplectic perturbations 
(i.e. perturbations that conserve all the Casimirs) is a necessary and sufficient condition 
of linear stability, it is clear that the minimization of E with respect to all perturbations 
that conserve S and I (i.e., two particular Casimirs) gives a sufficient condition of linear 
stability. In turn, this means that the maximization of S at constant E and I gives the 
same sufficient condition. 

Specializing Eq. (1261) to our case, we immediately obtain our sufficient condition: 



for all 8f(8,p) that at first order give 5S = 51 = 0, or equivalently 5E = 51 = 0. 

Another way to see that Eq. (179]) gives a sufficient condition, if Eq. ( J7"6l) gives a 
necessary and sufficient condition, is the following. Considering that De = 0, we have 
that all functions of the type 5f = Da give 5E = 51 = at first order. Therefore, we 
arrive at the conclusion that the condition for the maximization of S at constant E and /, 
or for the minimization of E at constant S and J, is stronger than the condition for linear 
dynamical stability. Said differently, if inequality ( 1791) is satisfied for all perturbations 
5f(6,p) that conserve E and I at first order, it is a fortiori satisfied for all perturbations 
that conserve E, I and all the Casimirs at first order. However, the reciprocal is wrong. 
Therefore, Eq. ( 1791) gives only a sufficient condition of linear stability. 

At this stage, it is interesting to note some analogies with thermodynamics. In 
particular, the formal stability obtained by maximizing S at constant E and I can 




(79) 
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be interpreted as a "microcanonical" formal stability problem if we regard S as a 
"pseudo entropy". Less refined stability properties can be found by relaxing one or 
both constraints (see below). 

On the other hand, taking S as being the Boltzmann entropy, we note that 
thermodynamical stability (in the usual sense) implies Vlasov linear dynamical stability. 
However, the converse may not be true in the general case, i.e. the Maxwell-Boltzmann 
distribution could be linearly stable according to f!76|) without being a maximum of 
Boltzmann entropy at fixed energy and normalization (i.e. a thermodynamical state). 

6.2. The "canonical" formal stability 

Following the usual procedure of thermodynamics, we pass from the "microcanonical" 
problem of maximizing S, Eq. (Q, at constant E and /, to the "canonical" problem 
of maximizing the "pseudo free energy" S — (3E (equivalent to minimizing E — ^S) 
at constant /. Introducing the Lagrange multiplier /i, we obtain again the first order 
variational problem: 

5S - (35E -n5I = Q, (80) 

equivalent to Eq. (fi~2l) . and therefore the same extremizing stationary state. Without 
repeating again the computations made in section [31 it is now clear that the condition 
of maximum (i.e., of formal stability) is given by the relation 

/oo r2n -i 

dp / dfl— — (5/(M) 2 (81) 
Jo 7{0,P) 

/oo / > 2tt poo r2n 

dp d9 dp' d9'6f(9 lP )V(9-9')5f(9',p')>0 
-oo JO J -oo Jo 

for all Sf(6,p) that at first order give 51 = 0. 

6.3. The "grand- canonical" formal stability 

Relaxing also the constraint of normalization is associated to the passage to the "grand- 
canonical" problem. Namely, we look for the maximum of the "pseudo grand-potential" 
S — f3E — jil (or the minimum of E — + ^1) without any constraint. The first order 
variational problem will be again given by Eq. flHUj) . thus obtaining the same stationary 
state, while the condition of formal stability is given by the relation 

/oo p2n -i 

dp / d9——(6f(9,p)) 2 (82) 
Jo i(Q,p) 

/oo p27t poo p2ir 

dp d9 dp' d9'5f(9,p)V(9-9')5f(9',p')>0 
-oo Jo J -oo Jo 

for all 5f(9, p). 

This unconstrained problem corresponds to the usual energy- Casimir method |28j . 
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6.4- Summary of stability problems 

We have found that the necessary and sufficient condition of linear dynamical stability 
for a stationary state f(e(9,p)) of the Vlasov equation is given by Eq. ( )66|) . We have 
proven that this is equivalent to the fact that the stationary distribution function / 
satisfies (locally) the problem: 

mm {E[f] | all Casimirs} . (83) 

This is the most refined criterion of formal stability as it takes into account an infinity 
of constraints. By relaxing some constraints, we have then found that progressively 
less refined, sufficient conditions of linear stability are given by the following problems. 
First, the "microcanonical" stability problem: 



max {S | E,I} , (84) 



equivalent to: 



min {E | S, 1} . (85) 

Then, the "canonical" stability problem: 

max{S-fiE | 1} . (86) 

Finally, the "grand-canonical" stability problem: 

max {S - @E - /il} . (87) 

The solution of an optimization problem is always solution of a more constrained 
dual problem [29]. Therefore, a distribution function that satisfies the "grand- 
canonical" stability problem (no constraint) will satisfy the "canonical" stability 
problem, a distribution that satisfies the "canonical" stability problem (one constraint) 
will satisfy the "microcanonical" stability problem, and a distribution that satisfies 
the "microcanonical" stability problem (two constraints) will satisfy the infinitely 
constrained stability problem flHBl . This is the analogous of what happens in the study 
of the stability of macrostates in thermodynamics. Of course, the converse of these 
statements is wrong and this is similar to the notion of ensembles inequivalence in 
thermodynamics. We have the chain of implications 

dHZD ^ dH6D ^ (JHSD ^ (JUD ^ flu. (88) 

The usefulness of these less refined conditions of linear stability will be clear after the 
stability conditions, that now appear as conditions to be satisfied by the perturbation to 
the stationary distribution function, will be transformed, in the application to the HMF 
model, in explicit conditions on the stationary distribution function itself. It will be 
shown that the less refined conditions of stability are associated to simpler expressions, 
and therefore, in a concrete calculation, one might use the simpler expressions if the 
more refined ones appear to be practically unfeasible. The procedure is to start by the 
simplest problem and progressively consider more and more refined stability problems 
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so as to prove (if necessary) the stability of a larger and larger class of distributions. Of 
course, if we can prove the stability of all the distribution functions with a particular 
criterion (see, e.g., Appendix), it is not necessary to consider more refined criteria. 

Remark: The connection between the optimization problems (!83|) - (|87|) was first 
discussed in relation to the Vlasov equation in [31], in Sec. 8.4 of [26] . and in Sec. 3.1 of 
[32] . Similar results are obtained in 2D fluid mechanics for the Euler-Poisson system [33J. 
Criterion fl83|) is equivalent to the so-called Kelvin-Arnol'd energy principle, criterion 
( |87l) is equivalent to the standard Casimir-energy method introduced by Arnol'd [M] 
and criterion (184"|) is equivalent to the refined stability criterion given by Ellis et al. [35J. 



7. The linear dynamical stability of Vlasov stationary states of the HMF 
model 



For the HMF model we have V(6 — 0') = — cos(0 — 0'). The extremization of a functional 
of the type ([9]) leads to a function of the type (see Eq. (fl6l)): 

n2 







p 



M x (f)cos6-M y (f)sm6\ + fi 



2 

with > 0, and with M x (f) and M y (f) given by self-consistency equations 



and 



M x {f) 



My(f) 



2tt 



dp' d6' cos 0'/(0', p' 



/■27T 

dp' d6' sin 07(0', p') . 
./o 



59) 



(90) 



(91) 



These are the two components of the magnetization. In this case the mean field potential 

mf) is 

$(0; /) = —M x (f) cos0 - M y (f) sin0 , 
and therefore the individual energy is given by: 



2 2 

e{e,p) = |- + $(0; /) = y - M x (f) cos0 - M tf (/) sin0 . 



(92) 
(93) 



Substituting Eqs. (192]). p)j) and (I9T1) in Eq. (J2SD, we obtain the total energy: 

^=y dp d0|-/(0,p) --(m, 2 (/)+m2(/)). 

The Vlasov equation for the HMF model reads: 
df(9,p) , 9/(0, p) 



+ p- 



-(M x (/)sin0-M 2/ (/)cos( 



9/(0, p) 







(94) 



(95) 



For the HMF model, we are interested in studying the stability of stationary solutions 
of Eq. f)5]) given by functions of the form /(0,p) = /(f - M x (f) cos0 - M B (/) sin0). 
Without loss of generality we can suppose that M y (f) = and therefore that 

,2 



/(M 



( y -Mcos0 j +// 



(96) 
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where for simplicity we have denoted M = M x (f) (dropping the explicit dependence on 
/). In this case we have: 

$(0;/) = -Mcosfl. (97) 

The individual energy is: 

2 2 

e(9,p) = |- + m f) = V - - Mcos0, (98) 



while the total energy is: 



poo p2n 2 i 

e = J dp de^f(e, P ) --m 2 



(99) 



The linearized Vlasov equation, governing the linear dynamics of 5f(6,p,t) around 
the stationary distribution f(6,p) of the form f[9"6l) . is obtained by linearizing Eq. (|95[) 
(in our case with M y (f) = and M x (f) = M), and it is given by: 

di 6f = ~ P W 6f + MSiT1 % 6f + Pf'«0, P ))^m 6 f) > ( 10Q ) 

where <&(9; Sf) must contain also the contribution of 5f to a magnetization in the y 
direction: 

/oo /*27r 
dp' d9' cos(6-6')5f(6',p'). (101) 
-oo JO 

Eq. f ll 001) is in the form of Eq. (13ip . with the linear operator D now taking the form: 

= PQg9(0,P) ~ M sm9— g{9,p) . (102) 

This operator has the property that De = and -D7 = 0, with e given in Eq. (1981) and 
where we again use for simplicity the notation r y{9 1 p) for the negative definite function 
f'(e). The linearized Vlasov equation can then be written, similarly to Eq. (13ip . as: 

= -j(9,p)(DK5f)(9,p,t) , (103) 
where for the HMF model the operator if is defined by: 

-1 /"OO /*27T 

(Kg)(9,p) = ——g(9,p)- dp' d9' cos(9 - 9')g(9',p') . (104) 

7(0, P) 7-00 Jo 

We are now in the position to follow the general results presented in section HJ 
Exploiting the particularly simple expression of the interaction potential in the HMF 
model we have, from Eq. (1661) . that the necessary and sufficient condition for the linear 
stability of f(9,p) is: 

d P r r d9-^—((D5f)(9,p)) 2 (105) 
Jo l{9,P) 



—00 

00 p2n \ z / poo p2n \ z 

dp J d9 cos 9(D5f)(9,p)\ -(J dp J d9 sin 9(D5f)(9,p)\ >0. 

We note that in Eqs. (1105P the first term in the left-hand side is positive definite, while 
the second and third terms are negative definite. 
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We have to find in which case the condition in Eq. (I105p is satisfied. We can 
exploit the antisymmetry of the operator D, that implies that the functional subspace 
orthogonal to the kernel of the operator is transformed in itself. In fact, if g\ belongs to 
the kernel of D, and g 2 is orthogonal to gi, then: 

(g 1 ,Dg 2 ) = -(Dg 1 ,g 2 )=0. (106) 

The kernel is made of the functions which depend on (9,p) through ^ — Mcos0^ . We 
may therefore transform the problem (11051) in the problem of satisfying the relation 

1 



dp 



2 77 



d9 



277 



7 (M 



(107) 



dp d9 cos 95 f (9, p) 



2 77 



dp d9 sin 05/(0, p) > 



subject to the conditions: 



J dp J d0^--Mcos0^ Sf(9,p) = 



0,1 



(106 



However, we should take into account the case in which the stationary distribution 
function / has a power law decay for large p. We therefore substitute the previous 
conditions with: 



dp 



277 



d9h 



V — - Mcos0 



y-Mcos0) Sf(9, p) =0 s 



0,1,... 



(109) 



where h is a function that assures integrability; it may be chosen, e.g., equal to 
exp - (£ -Mcos0) . 

Since the function 7(0, p) is even in 0, it is useful to separate 5f in its even and 
odd parts in 0, i.e. 8f(9,p) = Sf e (9,p) + Sf (9,p). In this way, our problem to satisfy 
the relation in Eq. (11071) subject to the conditions given in Eq. fll09j) . is separated in a 
pair of separate problems. Precisely, for the even part we have to satisfy: 



dp 



2- 



d0- 



-(5f e (9, P )) 2 - 



7(0, p) 

for all 5f e (9,p) such that the relations 



2k 



dp 



277 



d9h 



p 



— M cos 



p 

y 



dp d9 cos 95 f e {9,p) > (110) 



Mcos0 5f e (9,p) =0 s = 0,l, 



are verified. For the odd part we have to satisfy 



dp 



277 



d9- 



1 



-(Sf (9,p)f 



2- 



dp d9 sm95f (9,p) > (112) 



7(0, p) 
without any condition. 

At this point, we note the following things. Firstly, any 5f e such that cos05/ e has a 
vanishing integral will trivially give a positive value for the left-hand side of Eq. (IllOp . 
and similarly any 5f D such that sin05/ o has a vanishing integral will trivially give a 
positive value for the left-hand side of Eq. (I112p . Therefore the only problems can come 
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from function 8f e and Sf Q that do not have these mentioned properties. Secondly, since 
both (IllOp and (11121) are quadratic functions of 5f e and 5f Q , respectively, the sign of 
the expression is not changed by the multiplication of 5f e or 5f Q by any number. We 
can therefore study the sign of the left-hand sides of (11 101) and (II 12[) also by imposing 
a linear condition on 5f e and Sf Q . 

Therefore we proceed in the following way. We look for the extremum of the left- 
hand side of Eq. (IllOp . constrained by the conditions (II lip , with the further convenient 
constraint: 

r-2n 



dp d9 cos95f e (6,p) = 1. 



;ii3) 



o 



Similarly, we look for the extremum of the left-hand side of Eq. f l 1 X 2 j) under the 
constraint: 

r-2n 



dp d9 sm66f (6,p) = 1. 



114) 



o 



It is useful at this point to introduce the following definitions: 



a 



CO 



dp 



-oo 
oo 



2- 



d6h 



— -Mcos9 
2 



and 



dp 



dp 



2n 



d9h 



-oo 

OO 



2n 



d6h 



V 
~2 



- M cos 9 



Mcos9 



l(8,P) -Mcos9 
7 (M (e(0,p; f)) s 



(115) 



j(9,p) cos 9 



P 
~2 



Mcos9 



dp 



2jt 



d9h 



' 2 

— -Mcos9 



1 {9,p)cos9{e{9, F J)) 8 , (116) 



where the dependence on the function h is explicitly indicated. 

We begin with the problem related to 5f e . Introducing the Lagrange multipliers 
for the constraints (II lip and 2v for the constraint (1 1 1 3 j) . respectively, the 
conditioned extremum of the left-hand side of Eq. (11 10p is given by the equation: 



1 oo 

- ^rfpj 6 ^ P) - C 1 + v ) ^s9-J2 »i h) h 



s=0 



— - Mcos9 
2 



— -Mcos9 
2 



0.(117) 



It is clear that this extremum is a minimum, since the the second variation is simply 
— - > 0. Therefore the necessary and sufficient condition is that the disequality (IllOp 
is satisfied for the extremal Sf e (9,p) determined by Eq. (I117p . Denoting furthermore 
£ = -(l + i/), Eq. (HIID gives: 

oo 

Sfe(0,p) = £7(M cos^ - J2^h(0,p)h [e(9,p; /)] (e(9,p; f)) s . (118) 
Substituting this expression in Eqs. All ip . we obtain the system of equations: 

oo 

J2^' )a< s% = ^ h) 8 = 0,1,..., (119) 



s'=0 
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while substitution in Eq. ( II 13j) gives: 

/oo p2tt 00 

dp d9 1 (9,p)cos 2 9-J2^ ) vi h) -^- (120) 
-00 Jo s=0 

From the system d 1 1 9 j) . we may obtain the multipliers function of the multiplier 

£. We see in particular that the multipliers pO^ are proportional to £. We therefore 
introduce the "normalized" multipliers //£ , given by the solution of the system of 
equations: 

00 

£/# )a & = tf° * = 0,1,.... (121) 

s'=0 

We have that jjLp = ; substituting in Eq. dl 20[) we obtain: 

1 

dp d 7 (f, P~cos 2 5 - Er=o « 
The relation ( 11101) for 5/ e equal to the extremal function given by Eq. ( 11181) can 
now be easily obtained, taking into account Eqs. (11131) . f 1 1 2 1 j) and (I122p . Introducing 
the further short-hand notation 



oc 



J2ti h) 4 h) = zh) (123) 

s=0 

(where we put in evidence the dependence on the stationary distribution function 
through 7), we have: 

1 



.,_ > 1 • (124) 

z (i) - IZo d P Jo" d9 i( 9 iP) cos29 

We have thus obtained a relation involving only the stationary distribution function. 
This is the main expression of this paper. 

The relation valid in the case of the linear dynamical stability of homogeneous (i.e., 
with M = 0) stationary distribution functions is easily obtained. In fact, in that case 
r]i h ^ = 0; therefore jl^ = and thus 2(7) = 0. Then, taking into account that the 
integral of cos 2 9 is equal to 7T, Eq. fjl24j) becomes in this case: 



7T 



dp7(p)>0. (125) 



This is identical with the expression generally found in the literature for the linear 
stability of homogeneous distribution functions in the HMF model (see, e.g., Ref. 
[20| |22| |23j EH |26]); for this comparison we have to consider that for homogeneous 
distribution functions j(p) = f'(p)/p- We will show that for homogeneous distribution 
functions all extremal problems lead to the same result. This explains why the same 
expression is obtained from the "canonical" formal stability problem [23] (in principle, 
this approach only gives a sufficient condition of linear stability, but our present results 
show that it is in fact sufficient and necessary). 

The result just obtained for the problem associated to the even part Sf e (9,p) can 
be immediately transformed in that for the problem associated to the odd part Sf (9,p). 
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In this case, the only constraint is Eq. ( 11141) . so no multipliers fis are present. Then, 
we easily obtain the further condition, analogous to Eq. (I124p . that has to be satisfied 
by 7; namely: 

/OO /*27T 
dp / d9 -f(6,p) sin 2 9 > 0. (126) 
-00 Jo 

For homogeneous distribution functions this relation becomes equal to Eq. (I125p . 
However, it can be easily shown that for inhomogeneous distribution functions of the 
form given in Eq. ( 1961) . i.e., when M is strictly positive, Eq. H126I) is satisfied as an 
equality, independently of the particular form of the function and of the value of its 
parameters. In fact: 

/oo r2n roo i*2~k 

dp / d0 7 (0,p)sm 2 = dp d6 sin 2 6- 

■00 JO J -00 Jo 

1 pod p2n r\rp i poo i-2-k 

This equality is clearly associated to a Sf (9,p) that simply rotates, at first order, the 
distribution function f )96|) . This shows that, for inhomogeneous distribution functions 
(j96|), any odd 5f {6,p) will satisfy Eq. fTTT2]l . 

Summarizing the results of this section, the distribution function (I56"|) is linearly 
dynamically stable iff the relation (I124p is satisfied. This relation reduces to the simpler 
form given in Eq. (I125p for a homogeneous (i.e., with M = 0) distribution function. 

As we have shown in the general case, the most refined formal stability criterion 
( 1831) leads to the same necessary and sufficient condition. 



8. The formal stability of Vlasov stationary states of the HMF model: 
sufficient conditions of stability 

8.1. The "microcanonical" formal stability for the HMF model 

The problem related to the "microcanonical" formal stability is obtained by specializing 
to the HMF potential the expressions given in section 16.11 We thus obtain: 

dp 1^ d6— 1 (6f(9,p)) 2 (128) 
Jo 7(0, P) 

dp / d0 cos 05/(0, p)) - ( / dp d0 sin 05/(0, p) j >0 

-00 Jo J \J~oc Jo J 

for all 5/(0, p) such that the constraints of normalization and of total energy are satisfied 
at first order. Using the expression of the total energy for the HMF model, we have 
that the allowed Sf(9,p) have to satisfy: 

5£ = y dp jT d0 - M cos0^ 6f{9,p) = 0, (129) 

and 

POO f*27T 

51= dp d05/(0,p) = O. (130) 
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We immediately see that this is analogous to the linear stability problem, with the 
difference that now we have only the constraints associated to s = and s — 1, and the 
function h is the constant unitary function. In fact, the constraints (I130p and (I129p are 
nothing more than the constraints fillip for s = and s = 1, respectively, and in the 
case h = 1. Therefore, the only thing we need before writing down the result for this 
case is to adapt the definition of the parameters a s and rj s , and of the multipliers fl s , to 
the present situation. We thus define 

/OO /•27T 
dp I d9 1 (9,p)(e(9,p;f)) s (131) 
-oo 

and 



o 



/OO /•27T 
dp / d9 1 (9,p)cos9(e(9,p;f)) s . (132) 
-oo J0 

We therefore introduce the "normalized" multipliers fi S) given by the solution of the 
system of equations: 
i 

^2 ^s'Us+s' = rj s s = 0, 1. (133) 

s'=0 

The condition on / ~f(9,p) analogous to Eq. (I124p can now be immediately written 
down. It is given by: 

1 



w d) - IZc d p jlT^^'P) cos2 ^ 



> 1 , (134) 



where now the short-hand notation w(^) stands for: 



i 



^2^sVs = w{i). (135) 

s=0 

Summarizing, the distribution function (1§6"|) is formally stable with respect to the 
"microcanonical" criterion iff the relation (11341) is satisfied. This relation reduces to 
the simpler form given in Eq. (11251) for a homogeneous (i.e., with M = 0) distribution 
function, since in that case jl s = 0. Thus, dynamical linear stability and formal stability 
lead to identical conditions for stationary homogeneous distribution functions. 

Although it is not evident from the two expressions (11341) and (11241) . we know that 
if the former relation is satisfied, so is the latter, since we had found that the necessary 
and sufficient condition for the formal stability is also a sufficient condition for the linear 
dynamical stability. 



8.2. The "canonical" formal stability for the HMF model 

At this point, it is straighforward to derive the relation for the less refined formal 
stability problems. In particular, the "canonical" formal stability condition is completely 
analogous to the "microcanonical" formal stability condition, with the difference that 
the allowed 5f(9,p) in Eq. (11281) have to satisfy only the normalization constraint, i.e., 
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Eq. f )130p . This corresponds to taking only the constraint associated to s = 0. In 
particular, the system H133[) reduces to the single equation: 

h®o = Vo ■ (136) 
Then, the condition on 7(6*, p) now becomes: 

^ >1, (137) 

a-/-°°oo d f /o ?d0>y{0,p)co#0 

or, more explicitly, 

1 , X 

— —a >1- (138) 

— re*, a r2„ a ; dp L a ^ (0 , p )cos 2 

Again, for M = this condition becomes identical to Eq. (I125p . since in that case 
770 = 0. 

For consistency, it is interesting to note that the equality in Eq. (I138p corresponds 
to the condition of marginal stability found by another method in the Appendix F of 
[261. 



8.3. The "grand- canonical" formal stability for the HMF model 

Finally, the "grand-canonical" formal stability condition has no constraint at all 
Therefore, the condition on 7(0, p) is simply: 

1 



> 1 . (139) 

Again, for M = this condition becomes identical to Eq. (1125p . 

We remark the following point. The problem associated to the "grand-canonical" 
formal stability, Eq. ( 1871) . does not constraint the value of the normalization /, that 
therefore can be different from 1. It is clear however, that if we want to study the "grand- 
canonical" formal stability of a distribution function which is an extremum also for the 
other stability problems, we have to restrict ourselves only to normalized distribution 
functions. 

In the Appendix, just as a useful exercise, we show that, as we should expect, 
the inhomogeneous Maxwell-Boltzmann distribution function is "canonically" formally 
stable, and therefore also "microcanonically" formally stable and linearly dynamically 
stable, but it is not "grand-canonically" stable. This shows the role of the constraints 
and the importance of considering sufficiently refined stability criteria. 



9. Discussion and conclusions 



In this paper we have derived a necessary and sufficient condition of linear stability for a 
stationary state of the Vlasov equation. This condition is expressed by Eq. (I124p . which 
is the core of the paper. Less refined conditions of formal stability, which are sufficient, 
although not necessary, for linear dynamical stability, have also been obtained. 
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We should make two remarks. Firstly, it is clear that the form of the HMF 
interaction potential has simplified the task, and that further computations should 



dynamical stability will always require a degree of approximation, since the infinite sum 
implicit in the system ( 11211) and in the definition of 2(7), Eq. (11231) . will have to be 
replaced by some finite representation. This has led us to treat also the less refined 
formal stability conditions. At the price to have only sufficiency, more manageable 
expressions are to be expected. 

We would like to conclude with some comments about the relevance of the Vlasov 
stable stationary states from the point of view of thermodynamics. 

If the system is initially in a state that is not a stationary state of the Vlasov 
equation, it can be argued that there will be a rather fast evolution until a stable 
stationary state is reached. However, some care must be exercised about the sense in 
which this statement has to be taken. Analogously to the Liouville theorem for the 
iV-body distribution function of Hamiltonian systems, the time evolution of the one- 
body distribution function as governed by the Vlasov equation is such that its phase 
levels are conserved (in fact, the Vlasov equation states exactly the equality to zero of 
the convective derivative of /). In particular, an initial two levels /, i.e. an / which 
is constant in a given region of the (0,p) plane and zero outside of this region, will be 
two levels for all the following evolution. How can we expect such a function to evolve 
towards a smooth stable stationary state characterized by a continuity of phase values? 
This can be realized only in a coarse-grained sense, when we study a sort of smeared one- 
body distribution function, in which the value of / at each point is substituted by the 
average of / taken in a small neighbourhood of the given point. If there is an efficient 
mixing of the dynamics, we may expect that, no matter how small is the averaging 
neighbourhood (provided it is not vanishing), the averaged / will evolve towards the 
stable stationary state of the Vlasov equation. This is exactly the framework in which 
the Lynden-Bell theory of violent relaxation has been proposed [Ti] . 

It is a typical reasoning in thermodynamics or statistical mechanics to argue that 
the distribution functions of a system, in particular the one-body distribution function, 
will evolve according to the maximization of a functional given some constraints. For 
example, the final Boltzmann-Gibbs state of the one-body distribution function will be 
given by the maximization of the Boltzmann entropy 



subject to the constraints of normalization, Eq. (j3J), and given total energy, Eq. 
(JiU); the potential $(#; /) will have to be determined self-consistently. The use of 
the Boltzmann-Gibbs entropy (11401) is fully justified to characterize the state reached 
after the "collisional" regime has taken over; it will be the most mixed state given the 

* Interestingly, in stellar dynamics, using the Antonov criterion or the energy principle, it can be 



shown that all spherical galaxies with / = f(e) and /'(e) < are linearly [36], and even nonlinear ly 



be made for more 



complicated potentialo Secondly, actual computations for linear 





(140) 



[37] . stable. 
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constraints. In a collisionless regime such as the one governed by the Vlasov equation, 
one can make the same hypothesis about the evolution towards the most mixed state 
given the constraints, but one has to take into account that the mixing, even if maximally 
efficient, has to take place without violating the properties of the Vlasov equation, as the 
conservation of the phase levels of the distribution function. In this way, the Lynden-Bell 
expression of the entropy is obtained [T4"l [38]. 

The Lynden-Bell maximization problem gives a coarse-grained distribution of the 
form f(9,p) = f(e(9,p)) with / (e) < 0, i.e. a particular steady state of the Vlasov 
equation. As such, it extremizes a functional of the form 



at fixed normalization and energy. It can be shown that if this distribution function 
is a maximum of S at fixed I and E, then it is Lynden-Bell thermodynamically stable 
(see [33j EH] for the 2D Euler equation and Sec. V of |38j for the Vlasov equation). 
According to the present study, if it is a maximum of S at fixed I and E, it is also 
granted to be linearly dynamically Vlasov stable. More generally, it can be shown that 
the coarse-grained distribution function associated with a Lynden-Bell thermodynamical 
equilibrium is always dynamically stable (because it is a minimum of energy with respect 
to phase preserving perturbations), even if it is not a maximum of S at fixed E and I 
(see Sec. 7.8. of [33J for the 2D Euler equation). 

On the other hand, since it is not guaranteed that the mixing is always completely 
efficient (this is referred to as "incomplete relaxation"), one may argue that in cases 
when it is not efficient the dynamics will evolve trying to maximize, always in the coarse- 
grained sense, other functionals of the form (j!4ip that are not consistent with Lynden- 
Bell's theory (see [10] and Sec. XII of [38]). This is an essentially phenomenological 
approach. For example, the Tsallis functional is one particular case of such functionals 
(that are called generalized if -functions [40J ) . The constraints of normalization and 
total energy are always present, and it is immediate to see that, if the functional to 
maximize is of the form (j!4ip then the solution will always be a function of the form 

2 

f(9,p) = /(y + $(6*; /)). The problem at hand in this case is to show that a function of 
this form obtained by extremizing (j!4ip at fixed I and E is really a maximum, i.e., the 
value of (I14ip decreases if we perturb / without changing the values of the constraints. 
This "microcanonical" stability problem has been studied for the Tsallis distributions 



We are thus led to the conclusion that functions of the form f(9, p) = /(y + <I>(0; /)) 
are relevant also in a thermodynamical sense (with respect to the collisionless dynamics). 
However, we have shown that, while for homogeneous states formal "microcanonical" 
stability implies linear dynamical stability and viceversa, for inhomogenous states only 
the first implication is true. Thus, there might be inhomogeneous states which are 
formally "microcanonically" unstable (therefore not relevant from a thermodynamical 
point of view), that nevertheless are dynamically linearly stable. 




(141) 



in [27]. 



,2 



Stability of inhomogeneous quasi- stationary states 



30 



Appendix: The stability of the Maxwell-Boltzmann distribution function in 
the HMF model 

As an exercise we can apply the results of section [8] to the magnetized Maxwell- 
Boltzmann distribution function, that realizes the Boltzmann-Gibbs global equilibrium. 
For this case all the quantities can be obtained analytically, basically because in this 
case 7 = —/3f. We can consider from the beginning the "canonical" formal stability 
problem, given by Eq. (11371) . It is immediate to obtain in this case that ao = — and 
rjo = —f3M. We also have: 

/OO /*27T 
dp d 9 7(0, p) cos 2 9 = —f3 (M 2 + AM 2 ) , (142) 
-oo JO 

where we have denoted with AM 2 the variance of the magnetization, i.e., the expectation 
value ((cos 9 — M) 2 ). We then find that Eq. (11370 reduces to: 

/3AM 2 < 1 . (143) 

We first note that, in the homogeneous case, the last expression reduces to the known 
relation |/3 < 1, i.e., (3 < 2. However, it is not difficult to show that, when /3 > 2 and 
M > 0, Eq. (I143p is always satisfied, on the basis of the graphical construction that 
gives as a function of M, based on the relation 

where I\ and Iq are the modified Bessel functions of order 1 and 0, respectively. We 
then have: 

/3AM 2 = = _d_hWM]_ 1 

' P <9(/3M)/ (/3M) dMI (/3M) ' { ' 

where the last disequality is a consequence of the graphical solution of Eq. ( I144p . Then, 
we have proven the "canonical" formal stability, hence the "microcanonical" formal 
stability and the linear stability. On the contrary, it can be seen that the "grand- 
canonical" formal stability f 1 1 3 9 j) does not hold for magnetized states. In fact, from Eq. 
( I142p . we find that this stability would require 

/3M 2 + /3AM 2 < 1 , (146) 

that in turn, using Eqs. (I144p and (I145p . becomes: 



h(PM) \ 2 d h(/3M) 
I (/3M)J d(/3M)I (/3M) 



< 1 • (147) 



However, using the results in Appendix B of Ref. jH], it is proven that the left-hand 
side of the last expression, for M > is simply equal to /3 — 1. Since Boltzmann- 
Gibbs magnetized states are realized for /3 > 2, the disequality cannot be satisfied. If 
we consider a nonnormalized Maxwell-Boltzmann distribution function, and we denote 
with A its "mass", then the left-hand sides of the last two relations are multiplied by 
A. The stability condition then becomes a relation between A and /3: 

A(f3-l)<l. (148) 
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Finally, the reader might want to check that the "micro-canonical" formal stability 
condition (I134p . that must be satisfied since the "canonical" one is satisfied, reduces to 
the expression: 

1 - /3AM 2 + 2/3 2 M 2 AM 2 > . (149) 
Since Eq. (I143P is verified, then Eq. f!149j) is a fortiori verified, as it should. 
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